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Abstract 

SFOLD (Sfermion Full One Loop Decays) is a Fortran program package 
for calculating all sfermion two-body decay widths and the corresponding 
branching ratios at full one-loop level within the MSSM. The package adopts 
the SUSY Parameter Analysis convention and supports the SUSY Les Houches 
Accord input and output format. With the SFOLD package we found non- 
negligible electroweak corrections in bosonic decays of b, t and f. 
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External routines /libraries: SLHALib 2.2 
Subprograms used: LoopTools 2.6 
Nature of problem: 

If the MSSM is realized in nature, LHC will produce supersymmetric particles 
copiously. The best environment for a precise determination of the model param- 
eters would be a high energy e^e~ linear collider. Experimental accuracies are 
expected at the per-cent down to the per-mill level. These must be matched from 
the theoretical side. Therefore loop calculations are mandatory. 
Solution method: 

This program package calculates all sfermion two-body decay widths and the cor- 
responding branching ratios at full one-loop level within the MSSM. The renor- 
malization is done in the DR scheme following the SUSY Parameter Analysis 
convention. The program supports the SUSY Les Houches Accord input and out- 
put format. 
Restrictions: 

Unusual features: 

Additional comments: 

Running time: 

1. Introduction 

Supersymmetry (SUSY) is widely regarded as the most appealing ex- 
tension of the Standard Model (SM). Among supersymmetric theories, the 
Minimal Supersymmetric Standard Model (MSSM) is the most extensively 
studied. If the MSSM is realized in nature, LHC will produce supersymmet- 
ric particles copiously. The best environment for a precise determination of 
the model parameters would be a high energy e~^e~ linear collider. Exper- 
imental accuracies are expected at the per-cent down to the per-mill level 
[H, [2I, [sj . These must be matched from the theoretical side. Therefore loop 
calculations are mandatory. 

There are a few program packages available for the automatic computa- 
tion of amplitudes at full one-loop level in the MSSM: FeynArts/FormCalc 
H, SloopS 0,3 and GRACE/SUSY-loop [t"]. SloopS and GRACE/SUSY- 
loop also perform renormalization at one-loop level. However, so far there 
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is no publicly available code for the two-body sfermion decays at full one- 
loop level within the MSSM. Therefore, we have developed the Fortran code 
SFOLD [sj (and HFOLD Like HFOLD, it adopts the renormalization 

prescription of the SUSY Parameter Analysis project (SPA) [ll| and supports 



the SUSY Les Houches Accord (SLHA) input and output format p^. The 
package SFOLD (Sfermion Full One-Loop Decays) computes all two-body 
decay widths and the corresponding branching ratios of all sfermions at full 
one-loop level. 

Full one-loop radiative corrections to decays of sfermions into charginos 



and neutralinos are discussed in [13| for all sfermion flavours and generations. 
Yukawa corrections to sbottom decay into lighter stop and charged Higgs 
boson are given in [14|. SUSY-QCD corrections to top and bottom squark 



decays into all Higgs bosons are calculated in [15|. SUSY-QCD corrections 
to stop and sbottom decays into weak bosons can be found in 16|]. Finally, 
SUSY-QCD corrections to squark decays to gluinos are given in [l7|. Up 
to now, the electroweak corrections to sfermion decays into Higgs and gauge 
bosons have not been fully addressed. It turns out that also these corrections 
cannot be neglected in a significant part of the parameter space. 

2. MSSM sfermion sector at tree-level 
2.1. Masses and mixing angles 



The mass matrix in the gauge eigenstates basis {fi, /ij) is of the form [18 



m^? ajnif 



M)= ^7 (1) 



where 



= M?^ + (/f -e/s^)cos2^M| + mJ (2) 



7l Q,l 



7. 



R 



M^^^j, + efSyyCos2l3mz + mf (3) 



af = -/i*(tan^)-2^/^ (4) 

The mass eigenstates are obtained by diagonalizing A^^ with a unitary matrix 
= (cos 9j:, sin 9j; — sin Oj^, cos 9^) 

diag(mj-^ ,m^j^J = Rf M){R^y (5) 
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leading to the following sfermion masses and the mixing angle 
1 



m 



h 



2 V /i 



-afiTif 



ml — m?. Y 

' JL fl ' 



(6) 
(7) 



In the case of stop, sbottom and stau the left and right states are generally 
mixed. In contrast, sfermions from first and second generation have negligible 
Yukawa couplings. Therefore, /i = /l, /2 = Jr if < (A4p22 and 

/i = /r, h = h if (A^Pii > (A^P22. 

2.2. Decay patterns 

There are four possibilities of Feynman graphs for a two-body decay of a 
scalar: the decay into two scalars, into two fermions, into scalar and a vector 
particle and into two vector particles. The fourth possibility is not realized 
in the decay of a sfermion in the MSSM. The following sfermion decays are 
calculated (the first generation is shown, i,j,c— 1,2; n — 1, . . . , 4): 
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If the squark decay into a gluino is kinematically allowed it will dominate 
due to the QCD interaction. The third generation /2 can decay into fi 
and a neutral boson if there is sufficiently large mass splitting. For stops 
and sbottoms with large mass differences, decays into charged boson and a 
sfermion are possible. 



3. Calculation at full one-loop level 

The calculation of the decay widths is done in the same way as in the 
HFOLD program. We work in the DR (dimensional reduction) renormaliza- 
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tion scheme and in the general hnear gauge for the W and Z°-boson. 
We wrote a Mathematica program that generated the whole Fortran code 
using the packages FeynArts (FA) and FormCalc (FC). In FA all particle 
couplings in the MSSM are implemented. We kept the divergent parts of the 
counterterms to examine the UV finiteness of the renormalized amplitudes. 
The IR divergence is removed by using soft bremsstrahlung or by adding a 
corresponding 1 to 3 process (hard bremsstrahlung) for which we calculated 
all formulae analytically. The IR finiteness can be checked by varying the 
photon (gluon) mass A. Finally, we implemented several switchers that are 
described in sections 5.5, 5.6. 

4. Input parameters 

At the program start, SFOLD reads the file in SLHA format, where the 
Yukawa couplings, the gauge couplings gi, 92,93, gaugino masses, the soft 
breaking terms, the VEV, m^o, tan /3, fi are taken as input parameters at the 
scale Q. These parameters may be further changed. In that case, SFOLD 
recalculates on-shell masses of Susy particles and does not take them from 
the input file. 

5. Program manual 

5.1. Requirements 

• Fortran 77 (g77, ifort, gfortran) 

• C compiler (e.g. gcc) 

• LoopTools-2.6 [igl 

• SLHALib-2.2 [2Q] 

5.2. About version 1.0 

• The CKM matrix is set diagonal 

• Real SUSY input parameters 

• One-loop corrections to Higgs masses 

• Absence of three-particle sfermion decays 
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5.3. Installation 

1. Download the file sfold.tar.gz at 

|http: / / www.hephy.at / tools 

2. Unpack the archive by 

tar -xvzf sfold.tar.gz 

3. Go to the sfold folder and create symbolic links named LoopTools and 
SLHALib by 

In -s . .path. . /LoopTools-m.n LoopTools 
In -s . .path. ./SLHALib-m.n SLHALib 

4. Then run 

. / configure 
make 

5. That will generate an executable called sfold. To run SFOLD type 

./sfold 

5.4. Further notes 

• To use an older version of LoopTools you have to 

- change the 'call Itini' to 'call ffini' in decay. F file 

- include AOO, AOOC in looptools.h file 

• LoopTools and SLHALib are installed correctly if folders ix86-linux/bin 
and ix86-linux/include were created. (The name ix86-linux varies 
according to the system.) We further expect that the following files 
are present in the mentioned folders: looptools.h, libooptools . a, 
SLHA.h, libSLHA.a. 

5.5. The input file sfold. in 

1. name of input file (SLHA format) 

2. type = 1,2,3,4, 

1 = sneutrino, 2 = slepton, 3 = sup type, 4 = sdown type 

3. generation = 1,2,3 
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4. sfermion index = 1,2 

5. bremsstrahlung = 0,1,2 

= off, 1 = hard, 2 = soft 

6. resummation of bottom yukawa coupling = 0,1 

= off, 1 = on 

7. esoftmax 

cut on the soft photon (gluon) energy if soft strahlung is used 

8. name of output file 

5.6. The file sf old. F 

Further parameters, switchers and options are: 



delta_in 
lainbda_in 
Qscale 
xiW, xiZ 
local changes On 



osextmassesOn 
osloopmassesOn 



SMPRINTON 

MSSMPRINTON 

NO_SQUARK_MIXING 



divergent part of loop integrals 
photon (gluon) mass 

scale at which DR parameters are defined 
gauge parameters ^w,^z 

If set to 0, on-shell masses are taken from SLHA 
input file. If set to 1, on-shell masses are calculated 
through self energies. (Must be set to 1, if some of 
the input parameters described in section 4 are 
changed.) Masses can be seen in masses. out file. 
If set to 0, DR masses in kinematics. If set to 1, 
on-shell masses in kinematics. 



If set to 0, DR masses in vertex corrections. 

If set to 1, on-shell susy masses in vertex corrections. 

(This is to avoid a trap when a process is not 

kinematically allowed with respect to DR masses.) 

If defined, prints the SM parameters. 

If defined, prints the MSSM parameters. 

If defined, squark mixing is switched off. 



If any changes are done to sf old.F file, it is necessary to run make first. 
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5. 7. Working in Mathematica 

To work with SFOLD in Mathematica, the SPheno package (version 3.1.0) 
2l| is required additionally. To establish a link: 



1. Compile the sfold with Mmakefile by 

make -f Mmakefile (or make -f Mmakef ile-mac) 

2. Open the sf old.nb file 

3. Go to the terminal and type 

./Msfold -linkcreate & 

4. Copy the output, return to sf old.nb file and install link by 

link=Install [LinkConnect [" . . . output ..."]] 

5. Evaluate the cell. (It may be further necessary to go to the terminal 
and press Enter.) Continue working in the nb file. 

It is further neccesary to: 

1. Specify the path to the SPheno executable 

SPhenoPath = " . . .path. . ./SPheno" 

2. Setting the Directory to sfold directory 

SetDirectory[" . . .path. . ./sfold"] 

We provided the user with six examples. The first one calculates the partial 
widths of a sfermion at a specified mSugra point, the second one at a speci- 
fied MSSM point, the third one as a function of the mSugra parameter, the 
fourth one as a function of the GMSB parameter, the fifth one as a function 
of the AMSB parameter and the sixth one as a function of the MSSM pa- 
rameter. 

The input parameters are the same parameters as in sfold. in and sf old.F 
file except for the names of the SLHA input and output file which are set to 
be SPheno. spc, sfold. out, respectively. 

The following functions and auxiliary lists are implemented: 
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MakeSphenoSPCmsugra 
msugraPara 
MakeSphenoSPCgmsb 
gmsbPara 

MakeSphenoSPCamsb 
amsbPara 

MakeSphenoSPCmssm 
mssmPara 



Decay 
inputPara 



GetDM 



MakeSphenoSPCmsugra [s_String, l_List] 

{mO, ml2, tanb, sign(mu), AO} 
MakeSphenoSPCgmsb [s_String, l_List] 
{Lambda, Mines, tanb, sign(mu), n5, cgrav} 
MakeSphenoSPCamsb [s_String, l_List] 
{mO, m32, tanb, sign(mu)} 
MakeSphenoSPCmssm [s_String , l_List] 
{Q,Ml,M2,M3,Au,Ad,Ae,Ac,As,Amu,At,Ab, 
Atau,Mue,MA02,TB,MSLl,MSL2,MSL3,MSEl, 
MSE2,MSE3,MSQ1,MSQ2,MSQ3,MSU1,MSU2, 
MSU3,MSD1,MSD2,MSD3} 
Decay [LList] 

{type, gen, index,bremsOn,resumOn,esoftmax, 
Qscale,delta_in,lambda_in,xiZ,xiW,localchangesOn, 
osextmassesOn,osloopmassesOn} 
GetDM [l_List] 



MakeSphenoSPCmsugra [] , MakeSphenoSPCgmsb [] , MakeSphenoSPCamsb [] 
create the file SPheno.spc. At first, LesHouches.in file is created and then 
SPheno is called. s_String is the path to the SPheno executable, l_list is 
the list of the corresponding model parameters (see msugraPara, gmsbPara, 
amsbPara above). 

Decay [] calculates the partial widths of the specified sfermion. The rows 
correspond to decay modes, the first column is the tree level result, the 
second is SUSY-QCD and the third one is the full result. For the input 
parameter list, see inputPara above. 

GetDM [] returns the particular decay mode of the specified particle. l_List 
consists of the particle type, particle generation, particle index and of the 
decay mode index. 



6. Electroweak corrections in bosonic decays 

We give a few examples using the SFOLD package. First, we focus on 
the 62 decay. Figure [T] shows the partial widths as a function of the mSugra 
parameter mo. Other parameters are: nii = am^ + hrriQ + c, tan/3 = 3, 
sign(/i) = 1, = 0, where a, 6, c are chosen such that the whole parabola 
lies just above the excluded region in Figure 2 of reference 23|. The parabola 
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goes through the points [40, 330], [450, 300] and [740, 120]. For Figure [21 we 
took tan /3 = 10 and Aq = —300. Here the dominant decay is the decay into 
ii and W' if mo < 200. 
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Figure 1: S2 decays, tan/3 — 3,Ao ~ 0; red: txi, green: tx2J blue: 6x1, gold: tx2, 
violet: ixg, cyan: t^Sj, orange: bg, purple: ii W^. The solid lines correspond to the full 
one-loop result, dashed line to the SUSY-QCD result, dotted line to the tree level result. 



It is clearly seen that the electroweak corrections can reach about 20%. 
They cannot be neglected also in f2 Hq fi as follows from Figure El Here 
the f2 partial widths are functions of the MSSM parameter fi. Other MSSM 
parameters are: Mi = 100, M2 = 200, M3 = 600, = = 0, A^ = -400, 
A^ = As = 0, A^ = -400, At = -600, A = -900, A^ = -400, m^o = 140, 
tan/3 = 10, Ml = 250, = 100, Mq = Mj^ = M ^ = 500. 



7. Comparison of SFOLD with SPheno S.Obeta and SUSY-HIT 1.3 

In the following tables we compare the third generation sfermion partial 
decay widths with SPheno 3.0beta and SUSY-HIT 1.3 [22] at the SPSla' 
point. In SUSY-HIT, the QCD corrections to the decays g — + q', 
q ^ q' + and q ^ g + q are implemented. 
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r [GeV] 10 



200 



300 



400 




Figure 2: 62 decays, tan/3 — 10, = —300; red: txi, green: ^x^, blue: gold: tx2) 
violet: ixj], cyan: <X4, orange: 6g, purple: t\ . The solid lines correspond to the full 
one-loop result, dashed line to the SUSY-QCD result, dotted line to the tree level result. 




Figure 3: f2 decays; red: VtXij green: Pt X2 ^ blue: txi, gold: TX2, violet: TX3, cyan: 
TX4, orange: hofi, purple: fi Z. The solid lines correspond to the full one-loop result, 
dashed line to the SUSY-QCD result, dotted to the tree level result. 
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Or BR-tree SF-tree SF-sqcd SF-full SPheno SUSY-HIT 

v^xl 1-00 0.1166 0.1166 0.1124 0.1166 0.1099 

Table 1: Comparison of the partial decay widths of Ur 

fi BR-tree SF-tree SF-sqcd SF-full SPheno SUSY-HIT 

Txl 1-00 0.0166 0.0166 0.0161 0.0166 0.0123 

Table 2: Comparison of the partial decay widths of f\ 

fa BR-tree SF-tree SF-sqcd SF-fuII SPheno SUSY-HIT 

u^Xi 0.08 0.0148 0.0148 0.0151 0.0147 0.0089 

Txl 0.87 0.1548 0.1548 0.1481 0.1548 0.1513 

Txl 0.04 0.0080 0.0080 0.0080 0.0080 0.0044 

Table 3: Comparison of the partial decay widths of f2 

ti BR-tree SF-tree SF-sqcd SF-fuII SPheno SUSY-HIT 

txl 0.23 0.3023 0.3004 0.2901 0.3023 0.3139 

txl 0.05 0.0640 0.0674 0.0656 0.0640 0.0755 

hxt 0.72 0.9628 0.9747 0.9711 0.9771 0.1034 

Table 4: Comparison of the partial decay widths of ti 

h BR-tree SF-tree SF-sqcd SF-fuII SPheno SUSY-HIT 

txl 0.04 0.2518 0.2573 0.2235 0.2518 0.2664 

txl 0.10 0.6740 0.6197 0.6326 0.6740 0.6453 

txl 0.01 0.0732 0.0721 0.0733 0.0732 0.0886 

txl 0.04 0.2694 0.2818 0.2675 0.2694 0.3321 

bxt 0.26 1.7753 1.5696 1.6988 1.7485 1.6461 

bxt 0.15 1.0140 1.0280 0.9674 1.0223 1.0752 

h^h 0.05 0.3538 0.3729 0.2965 0.3528 0.4049 

tiZ 0.36 2.4851 2.5018 2.3873 2.4851 2.2822 

Table 5: Comparison of the partial decay widths of t2 
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0.36 


1.6161 


1.5735 


1.6665 


1.6078 


1.6682 


hxl 


0.04 


0.1621 


0.1550 


0.1319 


0.1621 


0.1610 


bxl 


0.29 


1.2901 


1.2086 


1.3346 


1.2901 


1.2810 


bxl 


0.00 


0.0104 


0.0103 


0.0111 


0.0104 


0.0112 


bxl 


0.00 
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Table 6 


Comparison of the partial decay widths of bi 
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0.1350 
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Table 7: Comparison of the partial decay widths of 62 
The screen output when running SFOLD is as follows: 



__ \l 1 
1 1 1 1 
1 1 1 1 


1 __ \ 
1 1 1 1 
1 1 1 1 


l__l 1 1 

/I 


_l l__l 1 

1 / 1.0 



Scalcir Full One Loop Decays by H. Hlucha, 
H. Eberl, W. Frisch 



FF 2.0, a package to evaluate one-loop integrals 
written by G. J. van Oldenborgh, NIKHEF-H, Amsterdam 



for the algorithms used see preprint NIKHEF-H 89/17, 
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'New Algorithms for One-loop Integrals', by G.J. van 
Oldenborgh and J. A.M. Vermaseren, published in 
Zeitschrift fuer Physik C46 (1990) 425. 



ffxdbO: IR divergent BO', using cutoff 1.00000000000000 
ffxcOi: infra-red divergent threepoint function, working with a cutoff 
1.00000000000000 

Qscale = O.lOOE+004 
lambda = O.lOOE+001 
delta = O.OOOE+000 

xiW = O.lOOE+001 

xiZ = O.lOOE+OOl 
DRbar parameters are read from SLHA input file SPheno-test . spc 
Masses are taken from SLHA input file. 
On-shell masses in kinematics. 
On-shell susy masses in vertex corrections. 
Hard photon (gluon) bremsstrahlung. 
BRs are written to SLHA output file outputs. slha 



tree 



~b. 
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-> 
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177440E+000 


/ BR 
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~b. 


.2 


-> 
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236841E+000 


/ BR 





21 


~b. 


.2 


-> 


b chi_20 





135022E+000 


/ BR 
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~b. 


.2 


-> 


b chi_30 





292376E-001 


/ BR 





03 


~b. 


.2 


-> 


b chi_40 





397420E-001 


/ BR 





04 


~b. 


.2 


-> 


~t_l w- 





514351E+000 


/ BR 





45 



Total width = . 113263E+001 



sqcd 



'b_2 -> t chi_l- 

'b_2 -> b chi.lO 

"b_2 -> b chi_20 

'b_2 -> b chi_30 



0.182192E+000 /BR : 0.16 

0.217873E+000 /BR : 0.19 

0.133641E+000 /BR : 0.11 

0.289576E-001 /BR : 0.02 
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~b_2 -> b chi_40 : . 395272E-001 / BR : 0.03 

~b_2 -> ~t_l W- : 0.562988E+000 /BR : 0.48 



Total width = 0.116518E+001 



full 



~b_2 


-> 


t chi_l- 





103772E+000 


/ 


BR 





13 


~b_2 


-> 


b chi.lO 





222152E+000 


/ 


BR 





27 


~b_2 


-> 


b chi_20 





806297E-001 


/ 


BR 





10 


~b_2 


-> 


b chi_30 





301927E-001 


/ 


BR 





04 


~b_2 


-> 


b chi_40 





365692E-001 


/ 


BR 





04 


~b_2 


-> 


~t_l w- 





340364E+000 


/ 


BR 





42 



Total width = . 813681E+000 
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